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STUDY QUESTION: Is there a contribution of the minor allele at the KRAS single nucleotide polymorphism (SNP) rs6 1 764370 in the /et-7 
microRNA-binding site to endometriosis risk? 

SUMMARY ANSWER! We found no evidence for association between endometriosis risk and rs6 1 764370 or any other SNPs in KRAS. 

WHAT IS KNOWN ALREADY: The rs6 1 764370 SNP in the 3' untranslated region of the KRAS gene is predicted to disrupt a com- 
plementary binding site (LCS6) for the let-7 microRNA, and was recently reported to be at a high frequency (31%) in 132 women of 
varying ancestry with endometriosis compared with frequencies in a database of population controls (up to 7.6% depending on ancestry), 
suggesting a strong effect of this KRAS SNP in the aetiology of endometriosis. 

STUDY DESIGN, SIZE AND DURATION: This was a case-control study with a total of I I 206 subjects. The study was performed 
between February 2012 and July 2012. 

PARTICIPANTS/MATERIALS, SETTING AND METHODS: We first investigated a possible association between common 
markers in KRAS and endometriosis risk from our genome-wide association (GWA) data in 3194 surgically confirmed endometriosis 
cases and 7060 controls of European ancestry. Although rs6 1 764370 was not genotyped on the GWA arrays, five SNPs typed in the 
study were highly correlated with this variant. The rs6 1 764370 and two SNPs highly correlated with rs6 1 764370 were then genotyped 
in 933 endometriosis cases and 952 controls using the Sequenom MassARRAY platform. 

MAIN RESULTS AND THE ROLE OF CHANCE: There was no evidence for an association between rs6 1 764370 and endometriosis 
risk P = 0.41 I and odds ratio =1.10 (95% confidence intervals: 0.88- 1 .36). We also found no evidence for an association between the 
highly correlated SNP rs 1 73870 1 9 and endometriosis. Their minor allele frequencies in cases and controls were of 0.087-0.091 similar 
to the population frequency reported previously for this variant in controls. Analyses of endometriosis cases with revised American Fertility 
Society stage lll/IV disease also showed no evidence for an association between these SNPs and endometriosis risk. 

LIMITATIONS AND REASONS FOR CAUTION: The GWA and genotyped data sets were not independent since individuals and 
cases from some families overlap. Controls in our GWA study were not screened for endometriosis. 

WIDER IMPLICATIONS OF THE FINDINGS: The key SNP, rs6 1 764370, was genotyped in a subset of samples. Our results do not 
support the suggestion that carrying the minor allele at rs6 1 764370 contributes to a significant number of endometriosis cases and 
rs6 1 764370 is, therefore, unlikely to be a useful marker of endometriosis risk. 

STUDY FUNDING/COMPETING INTEREST(s): The research was funded by grants from the Australian National Health and 
Medical Research Council and Wellcome Trust. None of the authors has competing interests for the study. 
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Introduction 

Endometriosis is a common, chronic gynaecological disease character- 
ized by the development of functional endometrial tissue outside the 
uterus. The disease is an important contributor to pelvic pain and sub- 
fertility. It affects 7- 1 0% of women of reproductive age and up to 50% 
of women with infertility (Treloar et al., 1 999; Giudice, 20 1 0). Medical 
treatment options are limited and surgical intervention may not 
prevent recurrence (Rogers et al., 2009; Giudice, 2010). 

Genetic factors contribute to endometriosis and it is inherited as a 
complex genetic trait (Kennedy, 1999; Treloar et al., 1999; Stefanson 
et al., 2002; Simpson and Bischoff, 2003). Numerous candidate gene 
studies have reported association for markers in candidate genes 
and endometriosis susceptibility, but results have generally not been 
replicated in subsequent studies (Simpson and Bischoff, 2003; Bischoff 
and Simpson, 2004; Guo, 2005; Montgomery et al., 2008). This lack of 
success is likely due to a number of reasons, including differences in 
disease definitions, case/control population sampling issues and the 
lack of power in small-scale studies, which detect only variants with 
large effects and are prone to detection of false-positive results (Zon- 
dervan et al., 2002). 

Two genome-wide association (GWA) studies have reported signifi- 
cant associations with endometriosis risk in three genomic regions on 
chromosomes 9p2 1 , 7p 1 5.2 and I p36 (Uno et al., 20 1 0; Painter et al., 
201 I). The regions showed strong evidence for association with 
disease risk and were replicated in independent data sets, demonstrat- 
ing that the studies had the power to detect common variants contrib- 
uting to disease risk. However, effect size estimates were small [odds 
ratios (ORs) < 1 .22] and the proportion of variance explained by the 
key single nucleotide polymorphism (SNP) on chromosome 7pl5.2 
(rs 1 2700667) was 0.36, or 0.69% of the estimated 51% heritability 
of endometriosis (Painter et al., 201 I). 

Recently, the minor allele at a SNP (rs6 1 764370) within the 3' un- 
translated region (UTR) of the Kirsten rat sarcoma viral oncogene 
homolog (KRAS) gene was reported to be at a high frequency (3 1%) 
in women with endometriosis, much higher than frequencies reported 
in unselected population samples (Grechukhina et al., 2012). KRAS is a 
strong candidate gene for endometriosis. The activation of Kras 
induced peritoneal endometriosis in 47% of mice, and 100% of mice 
developed benign endometriosis-like lesions within the ovarian 
surface epithelium (Dinulescu et al., 2005). Endometriosis-like 
lesions also developed and survived over a prolonged period in wild- 
type mice after transplanting endometrial tissue from mice with the 
activated Kras mutation (Cheng et al., 201 I). The rs6 1 764370 SNP 
is predicted to disrupt a complementary binding site (LCS6) for the 
/et-7 microRNA (miRNA) (Johnson et al., 2005; Esquela-Kerscher 
and Slack, 2006). The minor allele at rs6 1 764370 was shown to 
alter mRNA and protein expression levels in cultured endometrial 
stromal cells (Grechukhina et al., 2012), and stromal cells from 
women carrying the variant allele also showed increased proliferation 
and invasion. The authors concluded that carrying the minor allele at 
the rs6 1 764370 variant might account for up to one-third of endomet- 
riosis cases (Grechukhina et al., 2012). 

Previous studies have examined an association between variants in 
KRAS and endometriosis risk. Mutational screening of KRAS identified 
KRAS mutations in endometrioid carcinoma, but not in endometriosis 
lesions (Amemiya et al., 2004). Genotyping 32 tagging SNPs across 



KRAS in 958 endometriosis cases and 959 controls showed no evi- 
dence for common variants contributing to disease risk in a study 
that had 80% power to detect a risk variant with a minor allele fre- 
quency (MAF) of 5% at a genotype relative risk of 1 .7 (Zhao et al., 
2006). We would have expected the association signal for 
rs6 1 764370 to have been detected by correlation with tagging SNPs 
(Zhao et al., 2006) or in the GWA studies (Uno et al., 2010; 
Painter et al., 201 I). 

We therefore re-examined evidence for an association between 
common markers in KRAS and endometriosis risk from our GWA 
data (Painter et al., 20 1 I ) and genotyped rs6 1 764370 in endometriosis 
cases and controls. 

Materials and Methods 

Study samples 

We examined data from our previous GWA study including 3194 endo- 
metriosis cases and 7060 controls, with 2270 endometriosis cases from 
Australia (QIMR, Australian data set) and 924 cases from the UK 
(Oxford, UK data set) (Painter et al., 20 1 I ). The QIMR samples were un- 
related cases (one individual per family) drawn from 3908 affected women 
from multiple case families and single cases previously recruited into our 
endometriosis study. Women included in the study completed question- 
naires, provided a blood sample and gave access to their medical 
records, allowing retrospective confirmation of diagnosis for affected indi- 
viduals (Treloar et al., 2005; Painter et al., 201 I). UK cases included unre- 
lated cases drawn from 245 families with >2 affected sisters with surgically 
confirmed endometriosis (Treloar et al., 2002), and 785 singleton cases. 
We obtained medical records for all the cases and diagnosis of endomet- 
riosis was confirmed from surgical records in 100% of cases in Australia 
and 97% of cases in the UK (Painter et al., 201 I). Disease severity was 
retrospectively assessed using the revised American Fertility Society 
(rAFS) classification system (American Fertility Society, 1 985); A total of 
n= 1686 (52.7%) were diagnosed with stage A disease (rAFS l/ll or 
some ovarian disease with few adhesions), n = I 364 (42.7%) had stage 
B (stage lll/IV disease) and n = 144 (4.6%) had an unknown stage. 
In both the Australian and UK data sets, affected women completed a self- 
report questionnaire including questions about fertility history and pain 
symptoms; however, as questions were phrased differently, sub- 
phenotype analysis was limited to data from the largest case data set 
(QIMR, 2078 cases). 

The controls consisted of 1870 individuals recruited by QIMR and a 
further 5194 individuals provided by the Wellcome Trust Case Control 
Consortium 2 (WTCCC2) (Painter et al., 201 I). Ethics approval for the 
studies was obtained from the QIMR Human Ethics Research Committee, 
the Australian Twin Registry and the Oxford regional multi-centre and 
local research ethics committees. All the participants gave informed 
consent prior to testing. 

GWA association data 

Samples from the GWA study were genotyped on lllumina genotyping 
platforms. We reviewed a total of 97 common SNPs (MAF range from 
2.4 to 47.7%) located between 1 50 kb upstream and I 50 kb downstream 
of KRAS in all the endometriosis cases and in stage lll/IV cases compared 
with the controls. The KRAS variant rs6 1 764370 was not present on the 
chips, so we searched for correlated SNPs using the SNP Annotation 
and Proxy Search (SNAP) software available through the Broad Institute 
(http://www.broadinstitute.org/mpg/snap/ldsearch.php). Linkage dis- 
equilibrium (LD) is the non-random association of alleles at different loci 
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and indicates to what extent the SNPs are correlated in a particular ethnic 
population. Using the 1000 Genome pilot I data generated on 179 indivi- 
duals including 60 of European origin ('CEU'), the SNAP program com- 
pared the genotypes of SNPs within this region and identified those 
present on lllumina Human 610 and I M chips highly correlated 
(in strong LD) with rs6 1 764370. We identified five SNPs typed that 
were strongly correlated with rs6 1 764370 (r > 0.5) (Table I) and 
reviewed the association results previously calculated in our GWA study 
(Painter e£ al., 201 I). The effect of SNP genotype on the distributions 
of women with stage B disease or who answered 'Yes' to questions on 
subfertility and pain were tested by Pearson's chi-squared test. 

Genotyping 

To confirm results from the GWA data, we genotyped rs6 1 764370 and 
two correlated SNPs in the same sample of 958 endometriosis cases 
and 959 unrelated controls examined in our previous KRAS study 
(Zhao et al., 2006) using the Sequenom MassARRAY Genomics Platform. 
SNPs were genotyped as part of multiplex assays designed using the 
Sequenom MassARRAY Assay Design software (version 3.0: Sequenom 
Inc.) and samples were genotyped using standard methods (Zhao et al., 
2007, 2008). Departures from Hardy-Weinberg equilibrium and statistical 
analysis for association between SNPs and endometriosis risk were tested 
using the PLINK (http://pngu.mgh.harvard.edu/~purcell/plink/). Quality 
control of case samples from the GWA data analyses had identified 
case families with varying degrees of non-Caucasian ancestry. These 
were removed from the analyses of this data set and genotype data 
from 933 cases and 952 controls were included in the association analyses. 

Results 

We identified five SNPs present on lllumina HumanHap6 1 0-Quad and 
IM-Duo Chips that were highly correlated with the SNP rs6 1 764370 
genotyped in the Grechukhina et al. (20 1 2) study. These SNPs were in 
strong LD with rs6 1 764370 (r 2 > 0.568, Table I) and were all located 



within 140 kb of the 3'UTR for KRAS. One SNP rs 1 73870 1 9 was 
reported to be in perfect (r 2 = I ) LD with rs6 1 764370. There was 
no evidence for an association between this SNP and endometriosis 
risk (Table I) with P = 0.305 [OR = 0.95; 95% confidence intervals 
(CIs): 0.85-1.05]. The MAF for this SNP in our GWA data was 
0.087 in the cases, consistent with published estimates in other Cau- 
casian populations (Hapmap-CEU and EUR — interim phase I -1000 
Genomes) (Chin et al., 2008). 

There was no evidence for any association between endometriosis 
risk and other correlated SNPs, including rs859 14 1 , which was in high 
LD with rs6 1 764370 (r 2 = 0.789) (Table I) with P = 0.067 and OR = 
0.90 (95% CIs: 0.82- 1 .02). The latter SNP also has an MAF of 0.090. 
In addition, we checked all other common SNPs genotyped across the 
region, from 1 50 kb upstream to 1 50 kb downstream of KRAS in 
our GWAS data. There was no evidence for an association between 
any of the genotyped SNPs in the KRAS region and endometriosis 
risk (Fig. I). 

As cases in Grechukhina et o/.'s analysis included 89% with a rAFS 
stage lll/IV diagnosis, we restricted our analysis to cases diagnosed 
with rAFS stage lll/IV (stage B) disease, and the results showed that 
the minor allele at each correlated SNP was slightly lower in the 
cases compared with the controls (Table I). We also found no evi- 
dence for an association between endometriosis risk and any corre- 
lated SNPs (Table I). We further analysed the genotypic association 
between SNP rs 1 73870 1 9 and sub-phenotypes of endometriosis col- 
lected in the Australian data set (see methods). There was no evi- 
dence that genotypes at rs 1 73870 1 9 influenced the proportions of 
women with stage B disease, infertility or pain (Table II). 

To directly examine the association between rs6 1 764370 and 
endometriosis risk, we genotyped this SNP, together with 
rs 1 73870 1 9 and rs859l4l in our subset of endometriosis cases with 
a family history of disease and controls used in previous studies. 
There was no evidence for association between rs6 1 764370 and 



Table I Results of the tests of association for SNPs correlated with rs6 1 764370 in all endometriosis and stage B disease. 



SNP a 


Position 


r 2 with 

rs6 1764370* 


Al 


A2 


Allele frequency 
in cases 


Allele frequency 
in controls 


Chi-square 


P-value 


OR (95% CIs) 


All Cases 




















rs859l4l 


251 13291 


0.789 


G 


A 


0.090 


0.098 


3.346 


0.067 


0.90 (0.82- 1 .02) 


rs7303889 


25146242 


0.571 


C 


A 


0.173 


0.179 


0.877 


0.350 


0.96 (0.89- 1 .04) 


rs 173870 19 


25155526 


1.000 


G 


A 


0.087 


0.091 


1.053 


0.305 


0.95 (0.85-1.05) 


rs 17388893 


25285132 


0.568 


A 


C 


0.060 


0.064 


0.790 


0.375 


0.95 (0.84- 1 .07) 


rs 17329975 


25290239 


0.568 


C 


T 


0.060 


0.064 


1.145 


0.285 


0.93 (0.83-1.06) 


Stage B Cases 




















rs859l4l 


251 13291 


0.789 


G 


A 


0.094 


0.098 


0.441 


0.506 


0.95 (0.82-1.10) 


rs7303889 


25146242 


0.571 


C 


A 


0.179 


0.179 


0.000 


0.984 


0.99 (0.90-1.1 1) 


rs 173870 19 


25155526 


1.000 


G 


A 


0.089 


0.091 


0.080 


0.777 


0.98 (0.85-1.13) 


rs 17388893 


25285132 


0.568 


A 


C 


0.062 


0.064 


0.067 


0.795 


0.98 (0.83-1.16) 


rs 17329975 


25290239 


0.568 


C 


T 


0.061 


0.064 


0.302 


0.583 


0.95 (0.80-1.13) 



r : A measure of LD or non-random association for the observed frequencies of alleles at two markers measured as the square of the correlation coefficient. 
LD was assessed using the SNAP proxy search program (http://www.broadinstitute.org/mpg/snap/ldsearch.php). 
Al, reference allele; A2, alternative allele. 

a SNPs present on lllumina HumanHap6 1 0-Quad BeadChip and lllumina Human I M-DuoChip that are highly correlated with SNP rs6 1 764370 in 1000 Genomes Pilot I data of 60 Utah 
residents with European ancestry (CEU) individuals. 
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Figure I Evidence for association (— loglO P-values, /-axis) between genotyped SNPs and endometriosis across the chromosome 1 2 region ~ 1 50 kb 
from the 3'UTR KRAS variant (rs6 1 764370). In this Locus Zoom association plot SNP rs 1 73870 1 9 is represented by a purple diamond (http://csg.sph 
.umich.edu/locuszoom/) and is perfectly correlated with the 3'UTR KRAS variant rs6 1 764370. The red arrow shows the position of the 3'UTR KRAS 
variant rs6 1 764370. Other SNPs are colour coded according to the strength of their correlation (LD) with rs 1 73870 1 9 (measured by r 2 ). 



Table II Genotype at rs 1 73870 1 9 in the Australian 
cases and the number (and proportions) with 
moderate/severe (stage B) disease or who answered 
'Yes' to questions on subfertility or pain. 





Genotype at rs 173870 1 9 


P-value 


AA 


AG 


GG 


Stage B 


7 (0.39) 


149 (0.42) 


754 (0.39) 


0.74 


Subfertility 3 


5 (0.29) 


1 29 (0.37) 


717 (0.38) 


0.64 


Menstrual pain b 


17 (1.00) 


332 (0.94) 


1729 (0.91) 


0.13 


Pelvic pain c 


1 6 (0.94) 


294 (0.83) 


1540 (0.82) 


0.40 


Emergency 


10 (0.58) 


1 56 (0.45) 


792 (0.42) 


0.30 


treatment d 










Dysperunia 6 


1 1 (0.69) 


276 (0.79) 


1433 (0.77) 


0.48 



The number of women who answered 'Yes' to the following questions: 
a Have you tried for I 2 months or more on any occasion to conceive without 
success? 

Have you EVER experienced severe menstrual pain? 
c Have you EVER experienced severe pelvic pain? 
d Have you ever had to seek emergency treatment because of pain? 
e Have you ever experienced pain during sexual intercourse? 



endometriosis risk (Table III) with P = 0.4 1 I and OR = 1. 10 (95% 
CIs: 0.88-1.36). The MAF for this SNP was 0.099 in endometriosis 
cases. In agreement with the GWA data, we did not find any evidence 



for association between endometriosis risk and rs 1 73870 1 9 and 
rs859l4l (Table III). These two SNPs were also in high LD with the 
KRAS SNP rs6 1 764370 in our data set (r 2 = 0.858 and 0.537, respect- 
ively), although the estimates were slightly lower than in the 1000 
Genome pilot I data in our larger sample. 



Discussion 

Ras proteins have an important role as binary molecular switches in 
controlling signal transduction pathways. In mouse models, the activa- 
tion of Kras triggers the development of tissue resembling endometri- 
osis and endometrioid ovarian cancer (Dinulescu et al., 2005; Cheng 
et al., 201 I). In women with endometriosis, the frequency of the 
minor allele at rs6 1 764370 in the 3' UTR of KRAS gene was reported 
to be 3 1 % (4 1 / 1 32 cases), a much higher frequency than expected in 
the general population (Grechukhina et al., 2012). The authors con- 
cluded that rs6 1 764370 is a marker of endometriosis risk, potentially 
contributing to nearly one-third of all endometriosis cases and that it 
provides a novel method for diagnosis (Grechukhina et al., 2012). 
However, the results from our analyses in two large data sets 
provide no support for an association between rs6 1 764370 and endo- 
metriosis risk. 

The KRAS polymorphism rs6 1 764370 was not present on the lllu- 
mina 610 chip used for our GWAS. We therefore examined LD 
between rs6 1 764370 and the genotyped SNPs. SNP rs 1 73870 1 9 
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Table III Results of the tests of association for rs6 1 764370 and SNPs correlated with rs6 1 764370 in multiplex genotyping 
of endometriosis cases and controls. 



SNP 


Position 


r 2 with 
rs6 1764370 


Al 


A2 


Allele frequency 
in cases (n = 933) 


Allele frequency in 
controls (n = 952) 


Chi-square 


P-value 


OR (95% CIs) 


rs6 1 764370 


25360224 




G 


T 


0.099 


0.091 


0.677 


0.41 1 


1.10 (0.88-1.36) 


rs 173870 19 


25155526 


0.858 


A 


G 


0.089 


0.082 


0.509 


0.476 


1 .09 (0.86- 1 .37) 


rs859l4l 


251 1329 


0.537 


T 


C 


0.097 


0.090 


0.493 


0.483 


1.08 (0.87-1.35) 



LD was assessed using the PLINK program which was run on our multiplex genotyping data set. 



which was typed on the chips is in strong LD with rs6 1 764370 and 
genotyping of this SNP provides almost complete information on 
the other. We also identified four other SNPs typed on the lllumina 
610 chip and in moderate to high LD (r 2 > 0.57) with the KRAS 
3'UTR SNP. There was no evidence for an association with endomet- 
riosis for any of these correlated SNPs; allele frequencies in the cases 
and the controls were similar in our GWA sample including 3194 
endometriosis cases. 

Since rs6 1 764370 was not on the lllumina chips, we genotyped this 
SNP in a sample of endometriosis cases with a family history of 
disease. The studies were not independent since individual cases 
and some families overlap between the two sample sets. There was 
no evidence for an association as the allele frequency for 
rs6 1 764370 in our genotyped sample was 0.099 for the cases and 
0.09 1 for the controls, similar to the population frequency reported 
previously for this variant (Grechukhina et al., 2012). The results 
also agree with those of our previous study, where we found no asso- 
ciation between common variants in KRAS and endometriosis risk 
from genotyping 32 tagging SNPs across the KRAS gene (Zhao et al., 
2006). Analysis in our GWA data for 97 SNPs located within the 
region from 1 50 kb upstream and downstream of KRAS gene 
showed no evidence for different allele frequencies between the 
cases and the controls and no association signal between any SNP 
in KRAS and endometriosis risk. 

Of 132 women with endometriosis typed for the KRAS SNP 
rs6 1 764370 by Grechukhina et al. (2012), 89% (n = 117) of the 
cases were assigned stage lll/IV disease. When we limited our analysis 
to cases with stage B (lll/IV) disease, our data set showed no evidence 
of any association between endometriosis risk and SNPs highly corre- 
lated with rs6 1 764370 [proxy SNP rs 1 73870 1 9: OR = 0.98 (95% CIs: 
0.85-1.13)]. MAFs in stage B cases were slightly lower in the cases 
compared with the controls, very similar to those in all endometriosis 
cases and our Australian and UK control population samples. Taking 
these results together, our data provide no support for the minor 
allele of rs6 1 764370 conferring risk of endometriosis. 

Our analysis of KRAS variants is based on a sample of 2270 endo- 
metriosis cases recruited in Australia and 924 cases recruited in the 
UK (Painter et al., 201 I). The diagnosis of endometriosis was con- 
firmed from surgical records in 100% of cases in Australia and 97% 
of cases in the UK (Painter et al., 201 I). The controls in our GWAS 
were not screened for endometriosis. However, Grechukhina et al. 
used data from public databases for control allele frequency informa- 
tion and it is also unlikely that these individuals were screened for 
disease. 



In summary, we find no evidence for an association between risk of 
endometriosis and variants in KRAS including the rs6 1 764370 SNP in 
the let-7 miRNA-binding site in the 3'UTR of KRAS. Our results do 
not address the physiological consequences of carrying the variant 
allele in the let-7 miRNA-binding site. However, our study provides 
no support for suggestions that carrying the minor allele of 
rs6 1 764370 contributes to a significant number of endometriosis 
cases and rs6 1 764370 is unlikely to be a useful marker of endometri- 
osis risk. 
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